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Abstract 

Background: Real-time quantitative PCR (qPCR) is a technique frequently used to measure changes in mRNA 
expression. To ensure validity of experimental findings, it is important to normalize the qPCR data to reference 
genes that are stable and unaffected by the experimental treatment to correct for variability among samples. Unlike 
in some models of neuropathic pain, reference genes for models of inflammatory injury have not been validated. 
This study examined four candidate reference genes in an effort to identify and validate optimal genes for 
normalization of transcriptional changes occurring in the dorsal horn of the spinal cord and the rostral 
ventromedial medulla (RVM) following intraplantar injection of complete Freund's adjuvant (CFA). 

Results: The expression of hypoxanthine phosphoribosyltransferase 1 (Hprtl), beta-actin (Actb), mitogen-activated 
protein kinase 6 (Mapk6), and beta-2-microglobulin (B2m) was quantified in the dorsal horn and RVM of rats four 
days or two weeks after intraplantar injection of CFA or saline. The range of expression levels among these four 
genes differed by as much as 16-fold within the dorsal horn and the RVM. All four of these reference genes were 
stably expressed in both tissues and did not differ between saline and CFA-treated animals. Analyses using the 
statistical algorithms in geNorm and NormFinder programs determined that Mapk6 was the most stable gene and 
recommended the combination of Mapk6 and Actb, or Mapk6 and Hprtl, in such experimental conditions. 

Conclusions: This study validated the four genes Hprtl, Actb, Mapk6 or B2m and showed that any one or combination 
of two of them are good reference genes for normalization of mRNA expression in qPCR experiments in the spinal 
cord and RVM in the CFA model of inflammatory injury. 

Keywords: Complete Freund's adjuvant, Hprtl, Mapk6, Actb, B2m, Pain, Nociception, Rostral ventromedial medulla, 
Dorsal horn, Housekeeping genes 



Introduction sensitivity to both noxious and non-noxious stimuli that 

A better understanding of the molecular changes that can persist long after the resolution of inflammation 
occur during the development and maintenance of pain [1-4]. The underlying events entail transcriptional or 
after inflammatory injury is fundamentally important for translational changes in numerous genes and proteins 
the development of effective therapies. Peripheral inflam- including those involved in neurotransmitter release, 
matory injury activates signaling cascades in the peripheral receptor function and trafficking, subcellular signaling 
and central nervous systems as well as in the immune sys- pathways, and regulation of ion channel expression and 
tern. These cascades lead to tissue swelling and increased activity [1-4]. Identifying the nature of these changes can 

provide insights into the mechanisms of inflammatory in- 
jury and may elucidate new approaches for treating in- 



* Correspondence: donna-hammond@uiowa.edu flammatory pain. 



+ Equal contributors 

'Departments of Anesthesia, University of Iowa, Iowa City, IA 52242, USA 



Changes in gene expression that occur with peripheral 



'Pharmacology, University of iowa, 451 Newton Road, lowa City, ia 52242, inflammatory injury can be detected and measured in a 

USA 



o 



© 2014 Walder et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BlOlVICCl Central Commons Attribution License (http://creativecommons.Org/licenses/by/4.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain 
Dedication waiver (httpy/creati vecommons.org/publ icdomain/zero/1.0/) applies to the data made available in this article, 
unless otherwise stated. 



Walder et al. Molecular Pain 2014, 10:55 
http://www.molecularpain.eom/content/10/1/55 



Page 2 of 9 



sensitive and specific way using quantitative real-time 
polymerase chain reaction (qPCR) assays. Quantitative PCR 
is a frequently used technique that can be easily adapted 
to measure mRNA levels for any target protein and is the 
method of choice for absolute or relative quantification of 
mRNA expression. However, critics of qPCR maintain that 
it is often inadequately standardized and frequently in- 
consistent [5]. Some technical challenges inherent to the 
technique include isolating high quality RNA, identifying 
an efficient reverse transcriptase (RT) enzyme to generate 
cDNA, designing efficient and specific primers to amplify 
the desired mRNAs, and normalizing results to adequately 
validated reference genes. In 2009, Bustin et al. published 
the MIQE Guidelines: Minimum Information for Publica- 
tion of Quantitative Real-Time PCR Experiments in an ef- 
fort to standardize the information needed to ensure the 
relevance, accuracy, correct interpretation, and repeatabil- 
ity of qPCR experiments [6]. 

An essential component of qPCR is normalization of 
the target mRNA to a reference gene of interest in the 
same sample to control for variability associated with tem- 
plate input (amount of starting material) as well as RT and 
qPCR efficiencies. The reference gene mRNA should be 
stably expressed, and its abundance should show a strong 
correlation with the total amount of mRNA present in the 
sample [6]. Control reference genes were initially termed 
"housekeeping genes" because the genes were historic- 
ally chosen from cellular maintenance proteins that are 
ubiquitously expressed and whose mRNA was generally 
thought to have uniformly unchanging expression in dif- 
ferent cells and under different conditions. More recent 
reports have documented that housekeeping gene expres- 
sion can vary substantially (e.g. in tumor environments 
and other pathophysiological states) and that the choice of 
reference gene can significantly impact the conclusions of 
a study [7-9]. The development of statistical algorithms 
(e.g. geNorm, BestKeeper or NormFinder [10-12]) to help 
researchers determine the most stable reference genes and 
most appropriate combinations of reference genes has ad- 
vanced the field. 

Many studies of the molecular mechanisms that under- 
lie the development and maintenance of pain after inflam- 
matory injury have used qPCR to examine changes in 
gene expression in the peripheral and central nervous 
systems. Yet, to our knowledge, rigorous validation of ap- 
propriate reference genes has not been undertaken. This 
study quantified the expression of four potential reference 
genes, hypoxanthine phosphoribosyltransferase 1 (Hprtl), 
beta-actin (Actb), mitogen- activated protein kinase 6 
(Mapl<6), and beta-2-microglobulin (B2m), in the dorsal 
horn of the spinal cord and rostral ventromedial medulla 
(RVM) of the rat. These regions play important roles in 
the development and maintenance of nociception after 
peripheral inflammatory injury [1,4,13,14]. 



Results 

Levels of Hprtl, Actb, Mapk6, and B2m mRNA in the 
dorsal horn and RVM of rats with persistent inflammatory 
injury 

Transcripts for Hprtl, Actb, Mapk6, and B2m were 
readily detected in the dorsal horn of the spinal cord 
and the RVM of rats at four days (Figure 1A,C) and two 
weeks (Figure 1B,D) after intraplantar injection of saline 
or CFA. The four candidate reference genes were abun- 
dantly expressed in both tissues, with Cq values ranging 
from 21.11 to 21.85 for Hprtl, 17.81 to 18.42 for Actb, 
20.56 to 21.54 for Mapk6, and 19.07 to 19.44 for B2m. 
The relative abundance of these mRNAs in both the 
dorsal horn and the RVM was Actb > B2m > Mapk6 > 
Hprtl (P < 0.001). 

The levels of transcripts for all four candidate refer- 
ence genes in CFA-treated rats were comparable to 
those in the corresponding saline-treated group in the 
dorsal horn (Figure 1A,B; P > 0.5 for all genes at both 
time points) and in the RVM (Figure 1C,D; P > 0.3 for 
all genes at both time points). The expression of each of 
the four reference genes in the dorsal horn at four days 
did not differ from that at two weeks in CFA- or saline- 
treated groups (P > 0.3, each gene). The expression of 
B2m in the RVM also did not differ at four days and two 
weeks in CFA- or saline-treated groups (P > 0.2). The ex- 
pression levels of Hprtl, Actb, and MapK6 in the RVM 
at four days were statistically different from correspond- 
ing values at two weeks in saline- and CFA-treated rats 
(P < 0.05, each gene). However, the magnitude of the dif- 
ference was less than 0.2 Cq and unlikely to be biologic- 
ally significant. 

geNorm analysis of reference gene stability 

In the geNorm program, an M value under 1.0 is indica- 
tive of stable expression of mRNA [15]. In the dorsal horn, 
all four genes had M values less than 0.4, demonstrating 
that all four have high expression stability (Figure 2A,B). 
Comparisons of pair-wise variation in the stability values 
[12] of the four genes in the dorsal horn indicated that 
Mapl<6 and Hprtl were the most stable transcripts at four 
days and that Mapk6 and Actb were the most stable at 
two weeks. When the time points were pooled in the 
interest of following expression changes over time, geN- 
orm analysis indicated that Actb and Mapk6 were the two 
optimal genes in the dorsal horn (Figure 2C). Similarly, in 
the RVM, all four genes exhibited high stability of expres- 
sion at four days and two weeks after injury (Figure 2D,E). 
Comparisons of pair-wise variation in the stability value 
indicated that at four days Actb and Hprtl transcripts 
were the most stable, whereas at two weeks Hprtl and 
Mapk6 were the most stable (Figure 2D,E). Pooling data 
for both time points indicated that Hprtl and Mapk6 were 
the two most stable genes in the RVM (Figure 2F). 
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Figure 1 Determination of transcript levels (Cq) of the candidate reference genes in the dorsal horn of the spinal cord (A, B) and RVM 

(C, D). Levels were assessed four days (A, C) or two weeks (B, D) after intraplantar injection of CFA or saline. Values are the mean ± SEM 
of determinations made in 6-7 rats. In the dorsal horn, no differences were observed in the transcripts levels between saline- and CFA-treated 
rats, or between four days and two weeks, for any reference gene. In the RVM, transcript levels for each reference gene did not differ between 
saline- and CFA-treated rats. Although the Cq values at two weeks were statistically different from values at four days, the magnitude of 
the difference was about 0.2 Cq and unlikely to be biologically significant. 



NormFinder analysis of reference gene stability 

Expression stability for the four candidate reference genes 
was also calculated using NormFinder (Figure 3). In the 
dorsal horn, NormFinder analysis indicated that Hprtl 
and Mapl<6 were the best genes for normalization at four 
days and two weeks, respectively (Figure 3A,B). When the 
time points were pooled, NormFinder identified Mapk6 as 
the most stable gene for normalizing the expression data 
in the dorsal horn (Figure 3C). In the RVM, NormFinder 
analysis indicated that Actb and Hprtl were the best genes 
for normalization at four days and two weeks, respectively 
(Figure 3D,E). When the data for both time points were 
pooled for analysis, NormFinder identified Mapl<6 as the 
best gene for normalization in the RVM (Figure 3F). With 
respect to the best combination of two genes to normalize 
gene expression across time points, NormFinder identified 
Actb and Mapk6 in both the dorsal horn and RVM. 

NormFinder also includes an analysis of the optimal 
number of genes for normalization. In both dorsal horn 
and RVM, NormFinder identified the combination of 
Mapl<6, Actb, and Hprtl as optimal. However, the addition 
of one or two more genes caused at best a 0.06 decrease in 
the accumulated standard deviation, much less than the 
0.15 that is considered a significant change (Figure 4). 
Thus, although the use of all three genes is predicted to be 



optimal, one gene would be acceptable for normalization 
of a target gene when resources are restricted. 

Discussion 

As qPCR methodology has matured and its use has be- 
come more widespread, an increasing number of reports 
have demonstrated that the expression of classically used 
reference genes can vary substantially depending on the 
experimental conditions [7-9,16]. Moreover, radically 
different conclusions can be reached depending on the 
reference gene used for normalization [16-19]. A PubMed 
literature search (March 20, 2014; search algorithm avail- 
able from authors) identified 260 papers in English that 
described qPCR in rat models of nociception. The major- 
ity of the papers (83%, n = 216) normalized the expression 
of the target genes to a single reference gene. Glyceral- 
dehyde 3-phosphate dehydrogenase (Gapdh) was used 
in 48% (n = 104), Actb in 25% (n = 54), Hprtl in 4% 
(n = 8), and other housekeeping genes in the remaining 
23% (n = 50) of the papers. Only 5% (n = 12) of the papers 
used more than one reference gene for normalization. 
When the search was restricted to studies using the CFA 
model of inflammatory injury in the rat, only 16 studies 
were identified. These 16 papers also normalized qPCR 
data with the same reference genes and in the same 
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Figure 2 Assessment of reference gene stability using the geNorm algorithm. Average expression stability (M) for genes in the dorsal horn 
of the spinal cord (A, B, C) and in the RVM (D, E, F), at four days (A, D) and two weeks (B,E) after intraplantar injection of CFA or saline, or for 
both time points combined (C, F). All four genes were very stably expressed (M < 1 .0) in both tissues and at all time points. The lowest M value 
indicates the most stable gene. 



relative ratios as observed in the total 260 papers. To our 
knowledge, the present study is the first to specifically val- 
idate the use of reference genes for normalizing qPCR 
data in the CFA model of inflammatory injury. 

Optimally, the transcript of the reference gene(s) se- 
lected for any given qPCR analysis will be expressed in an 
abundance comparable to that of the mRNA of interest 
[6]. This study identified four different reference genes - 
Mapk6, Actb, Hprtl, and B2m - that can be used for the 
normalization of qPCR data in the CFA model of inflam- 
matory injury in rats. These four genes were selected be- 
cause they are broadly expressed, involved in different 
homeostatic functions (Table 1) and are not co-regulated 
[12]. Moreover, two of these genes were also studied for 
validation as potential reference genes in the spared nerve 
injury [20] and spinal nerve ligation [21] models. The 
addition of our data allows the assessment of the stability 
of those genes over different pain models of chronic pain, 
possibly strengthening their utility as reference genes. The 
Cq values of the four genes ranged from 17 to 21.85 Cq, 
reflecting a 2 or a 16-fold difference in expression. Each 
gene was stably expressed within each time point: four 
days (subacute injury) and two weeks (chronic injury) after 
injury in CFA-treated animals compared to corresponding 
saline-treated animals. Interestingly, the Cq values for 



these genes in the RVM were very similar to those in the 
dorsal horn. The NormFinder algorithm indicated that 
Mapl<6 was the reference gene with the least variability 
and therefore likely the best reference gene of the four in 
our experimental conditions. It also indicated that the 
combined use of three reference genes (Mapk6, Actb and 
Hprtl) was optimal for normalization in both the dorsal 
horn and the RVM. However, the gain in stability between 
the use of a single gene and the use of three genes is min- 
imal, as shown by the very small changes in the accumu- 
lated SD (< 0.06) between one and three genes, indicating 
that any one of these three genes would be acceptable as a 
unique reference gene. 

The need to report qPCR data in a standardized format 
is now widely accepted, and recommendations for the uni- 
formity and reproducibility of qPCR experiments are listed 
in the MIQE Guidelines [6]. These guidelines emphasize 
the need to control for sample-to-sample variation by 
normalization with reference genes. Normalization is ne- 
cessary for reliable qPCR studies because the starting 
material, RNA extraction, RT efficiency, and qPCR effi- 
ciency can vary among experiments. Moreover, gene ex- 
pression is highly tissue-specific and often varies based on 
the pathophysiological status of the organism or experi- 
mental treatment. Thus, it is imperative that a pilot study 
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Figure 3 Assessment of reference gene stability using the NormFinder analysis. Stability values for genes in the dorsal horn (A, B, C) and 
in the RVM (D, E, F), at four days (A, D) and two weeks (B, E) after intraplantar injection of CFA or saline, or for both time points combined 
(C, F) are shown. The lowest standard deviation (SD) indicates the most stable gene. 



be conducted at the outset of any study to identify the 
optimal reference gene or combination of genes for that 
specific experiment. Unfortunately, very few validation 
studies for the use of reference genes exist in the pain 
field, and all are restricted to nerve injury models. Several 
studies have identified suitable qPCR reference genes 
in dorsal root ganglia in models of neuropathic pain 
[11,21,22]. For example, Mapk6 and Gapdh were identi- 
fied as the two most stably expressed genes in an analysis 
of L4 and L5 dorsal root ganglia following L5 spinal nerve 
ligation [21]. Hprtl and 18S were validated as stable genes 
for normalizing expression levels in the dorsal root ganglia 
in the spared nerve injury model [22,23], although this 
conclusion was not shared by Piller et al. [20] . Seven com- 
monly used reference genes, of which Actb was the most 
stable, were ranked and validated as good reference genes 
in spinal cord tissue in the spared nerve injury model of 
neuropathic pain [20]. These reports support the conclu- 
sion that the optimal reference gene(s) will be specific to 
each experimental condition, and likely gender as well. 

Gapdh was the single most frequently used reference 
gene in the literature search but was not included in this 
analysis for several reasons. Over the last decade, the 
transcription of Gapdh has been reported to be signifi- 
cantly regulated in different experimental settings and is 
variable in different tissues [24-26]. Moreover, the rat 
genome contains 329 Gapdh pseudogenes [27,28], some 



of which are transcribed and have the same sequence as 
the active Gapdh transcript and, as such, can be detected 
by primers and amplified. For example, amplification 
using a published Gapdh primer pair [29]: forward (5'- 
ACCACGAGAAATATGACAACTCCC) and reverse (5'- 
CCAAAGTTGTCATGGATGACC), designed to amplify 
a product of 100 base pairs, can give the same size 
amplicon from numerous transcripts. Performing an in 
silico PCR test by inserting the primer pair sequences 
for Gapdh in the UCSC genome website revealed 53 
PCR products, each with 100 base pairs and 18 add- 
itional matches ranging from 75 to 120 base pairs, de- 
rived from genes on 16 different chromosomes. Other 
published primers for Gapdh yielded similar findings 
[20,21]. Even if no genomic DNA contamination exists, 
and the RNA samples are treated with DNase, a more op- 
timal qPCR assay design is one that does not amplify any 
pseudogenes. Similarly, Actb has pseudogenes [28], but 
the primers used in the present study were carefully se- 
lected to amplify Actb cDNA and not any of the pseudo- 
gene sequences. The UCSC in silico PCR search identified 
only one unique product in the rat genome with the Actb 
primers used here. Hprtl also has pseudogenes, but the 
primers selected for this study are unique: one primer 
spans the splice junction of the gene, and when the primer 
pair is tested in silico, no genomic product was identified 
since the primers can only amplify the mRNA. 
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Figure 4 Determination of the optimal number of reference 
genes by the NormFinder algorithm. Stepwise calculation of 
stability began with the most stable gene (Mapk6), followed by 
serial addition of the next most stable gene and calculation of the 
accumulated standard deviation (circles). In both the dorsal horn 
(A) and in the RVM (B), the optimal number of reference genes to 
be used in qPCR experiments was three (open circle). 



Conclusions 

In summary, this study identified four reference genes that 
are stably expressed in the dorsal horn and the RVM four 
days and two weeks after intraplantar injection of CFA or 
saline. Each assay for the four reference genes was de- 
signed with primers that were unique for their target and 
would not amplify pseudogenes. Mapl<6 was the most 
stable reference gene, although all genes tested met the cri- 
teria for a validated reference gene. We suggest that using 
any one or combination of two of these assays (e.g. Mapk6 
and Actb, or Mapk6 and Hprtl) for normalization of data 
would yield accurate and reproducible results when study- 
ing mRNA expression in the spinal cord and RVM in the 
CFA model of inflammatory injury. 

Methods 

Animals and inflammatory injury model 

Male Sprague-Dawley rats (200-350 g; Charles River, 
Raleigh, NC, USA) were housed in pairs in the University 
of Iowa Animal Care Facility in rooms with a 12 hr light/ 
dark cycle with water and food provided ad libitum. All 
studies were conducted in accordance with the Guide for 
Care and Use of Laboratory Animals published by the 
National Institutes of Health following the guidelines 
of the International Association for the Study of Pain. 
The experiments were approved by the University of Iowa 
Animal Care and Use Committee (protocol 1107156), and 
care was taken to minimize the number of animals used 
and their suffering. 

Persistent inflammatory injury was induced by intra- 
plantar injection of CFA. Adult rats were lightly anesthe- 
tized with isoflurane, and the thickness of the hindpaw 
was measured with digital calipers. The plantar surface 
of the left hindpaw was injected with 150 ul of CFA 
(150 ug of Mycobacterium butyricum, Calbiochem, La Jolla, 
CA, USA) or sterile-filtered saline at pH 7.4. The rats were 
returned to their cages and singly housed for four days or 
two weeks, depending on the experiment. After the rats 
were euthanized, the thickness of the hindpaw was mea- 
sured to verily the presence of inflammation. Four days 
after injection of CFA, the thickness of the ipsilateral hind- 
paw had increased from 6.1 ± 0.1 to 9.4 ± 0.4 mm (N = 12; 



Table 1 Primer sequences for the candidate reference genes 



Gene 
name 


Accession number 


Function 


Sequence 


Product size 
(base pairs) 


Hprtl 


Genbank NM_012583 


Purine synthesis in salvage pathways 


Forward 5'-CTCATGGACTGATTATGGACAGGAC 
Reverse 5'-GCAGGTCAGGAAAGAAGTTATAGCC 


123 


Actb 


Genbank: NMJB1144 


Cytoskeletal structural protein 


Forward 5 '-CGGCG AGTACAACGTTCTTG 
Reverse 5'-GCAGGGATATCGTCATCCAT 


81 


Mapk6 


Genbank: NM_031622 


Member of the Ser^hr protein 
kinase superfamily 


Forward 5 '-TAAAGCCATTGACATGTGGG 
Reverse 5'-TCGTGGAGAACAGGGATAGA 


129 


B2m 


Genbank: NM_012512 


Beta-chain of major histocompatibility 
complex class I molecules 


Forward 5'-CGAGAGCGATGTATATGCTTGC 
Reverse 5'-GTCCAGATGATTCAGAGCTCCA 


114 
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P < 0.01). The thickness of the hindpaw in the saline-treated 
cohort was unchanged (6.2 ± 0.1 mm; N = 13; P > 0.3). Two 
weeks after injection of CFA, the thickness of the ipsilat- 
eral hindpaw had increased from 6.0 ± 0.1 to 9.0 ± 0.3 mm 
(N = 10, P < 0.01) whereas the thickness of the hindpaw 
in the saline-treated cohort was unchanged (6.0 ± 0.1 mm, 
N = 10, P > 0.4). Measures of nociception were not made 
to minimize stimulation. 

qPCR - Quantitative Real-Time PCR Analysis 
Tissue dissection 

On the designated day, the rats were euthanized by C0 2 
inhalation and a 2-mm transverse slice of the brainstem 
containing the RVM was immediately isolated on ice and 
frozen on a platform on dry ice. To obtain the RVM, a 
1.5-mm diameter tissue punch (Harris Unicore, Ted Pella 
Inc., Redding, CA), centered on the midline immediately 
above the pyramids, was removed from the frozen slice of 
brainstem tissue. The remainder of the slice was fixed in 
10% formalin containing 30% sucrose to allow the verifica- 
tion of the site of the RVM tissue punch. The L4 and L5 
portion of the spinal cord was removed, chilled on an ice- 
cold platform, and the ipsilateral dorsal horn was excised. 
The tissues were stored at -20°C in RNAlatef" (Ambion, 
Life Technologies, Carlsbad, CA) until RNA isolation. 

Primer design 

DNA oligonucleotide primers were synthesized and pur- 
chased from Integrated DNA Technologies (Coralville, IA). 
The sequence of the forward and reverse primers for 
each of the four reference genes are listed in Table 1. 
The primers for Mapk6, B2m, Hprtl, and Actb are 
described elsewhere [21,30,31], and checked with the 
Primer 3 software (http://biotools.umassmed.edu) [32]. 
Each qPCR assay consists of primers that hybridize to se- 
quences that lie on different exons or span a splice junc- 
tion, separated by one or more introns, such that when 
qPCR is conducted, only the cDNA sequence is ampli- 
fied. The Ensembl database (http://useast.ensembl.org/) 
was used to examine the genomic structure of the gene 
and its transcripts. Special care was used in the selection 
of primers for Actb and Hprtl, which are known to have 
pseudogenes [12,27,28]. Contaminating genomic DNA, 
if present, will not amplify any product under the re- 
action conditions. Thus, only the specific mRNA tar- 
gets will be measured. When the forward and reverse 
primer sequences are entered into the UCSC genome bio- 
informatics website (http://genome.ucsc.edu/) and tested 
in silico against rat genomic DNA, only a single product is 
identified for Mapk6 and Actb primers. In the case of the 
Hprtl and B2m primers, no genomic match is found be- 
cause one of primers spans a splice junction and can only 
hybridize with the correct cDNA sequence. In addition, 
the PCR amplicons for each qPCR assay were cloned into 



pSC-A, the PCR cloning vector, according to the manufac- 
turer's protocol (StrataClone PCR Cloning Kit, Agilent 
Technologies, Santa Clara, CA), and sequenced at the 
Iowa Institute of Human Genetics, Genomics Division. 
Results were aligned with the Genbank sequence for the 
intended mRNA using the NCBI BLAST program to con- 
firm the specificity of each primer pair. When the primers 
are used in a qPCR assay, they amplify a unique species 
with a sharp melting curve, and no primer-dimer products 
are detected. 

RNA isolation, RT, and qPCR 

Total RNA was isolated from dorsal horn and RVM tissue 
according to the manufacturer's protocol (RNeasy Lipid 
Tissue Mini Kit, Qiagen). Briefly, each tissue sample was 
homogenized in 1 ml of QIAzol lysis reagent. The lysate 
was extracted with chloroform, centrifuged, and the 
supernatant saved in a clean tube. The supernatant frac- 
tion was mixed with an equal volume of 70% ethanol 
and loaded on the column. Samples were treated with 
DNase I while on the column for 15 min at room 
temperature. After washing, the RNA was eluted from 
the column using RNase-free water. The concentration 
of total RNA was measured on a Nanodrop spectropho- 
tometer (ND1000 3.8.1, Thermo Scientific, Wilmington, 
DE). The RNA integrity number was determined for 20% 
of the RNA samples as further RNA quality control and 
routinely showed values > 9.0 (out of 10) (Agilent Model 
2100 Bioanalyzer, Santa Clara, CA). Reverse transcription 
was performed according to the manufacturer's protocol 
using 100-600 ng of purified RNA and the Superscript 
VILO cDNA synthesis kit (Life Technologies, Carlsbad, 
CA) in a 20 ul reaction volume. No-reverse transcript- 
ase controls were also run for each RNA sample. The 
qPCR was performed with each well containing the 
cDNA product generated from 5 ng of input RNA, for- 
ward and reverse primers (12 nM), and the fluorogenic 
DNA-binding dye iQ™ SYBR 8 Green Supermix (Bio-Rad, 
Hercules, CA) in 20 ul. Reactions were performed in 
triplicate on a Bio-Rad CFX96 thermocycler (Bio-Rad, 
Hercules, CA). The cycle conditions were: 50°C for 
2 min, 95°C for 10 min followed by 40 cycles of (95° 
C for 15 s, 60°C for 1 min and 72°C for 1 min), and 
then 95°C for 1 min and 55°C for 1 min. A thermal 
melting curve was generated from 55 to 95°C, at in- 
crements of 0.5°C for 10 s. No reverse transcriptase 
and no template controls for each primer pair were 
also tested in triplicate and did not amplify any prod- 
uct. Amplification efficiencies, calculated using the 
Bio-Rad CFX Manager 3.0 software, were similar for all 
primers and averaged at E = 101.2% ± 2.8, r 2 > 0.993 ± 
0.002, slope = -3.30 ± 0.07. Cq values for each sample 
were calculated by the Bio-Rad CFX96 software (Bio-Rad, 
Hercules, CA). 
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geNorm and NormFinder 

These two programs are available in the GenEX 6 soft- 
ware package (http://genex.gene-quantification.info/). The 
geNorm program calculates gene expression stability (M) 
for each reference gene as the average of a pairwise vari- 
ation for the reference gene relative to the others. The 
stepwise analysis allows for a ranking of the genes accord- 
ing to the calculated M value. The lower the M value, the 
more stable the reference gene. 

The data for NormFinder are organized in groups: 
e.g. A) four day RVM, B) two week RVM, C) four day 
dorsal horn, and D) two week dorsal horn. The program 
calculates intra-group and inter-group variations in the 
expression of the reference genes from which a stability 
value is generated for each reference gene. The candidate 
gene with the lowest stability value is considered the most 
stable reference gene. By calculation of the standard devi- 
ation of the data, the NormFinder program selects the 
best number of reference genes for an experiment. 

Statistical analysis 

Data are expressed as mean and S.E.M. The rank order of 
the expression of the four genes was analyzed by a one- 
way ANOVA followed by Bonferroni's multiple compari- 
son post-hoc test. A two-way ANOVA (factors: treatment 
and time) was used to compare the Cq values between 
saline- and CFA-treated rats, as well as between the four 
day and two week time points for each gene within each 
tissue. A P < 0.05 was considered significant. 
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